Topological Features in Ion Trap Holonomic Computation 
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Topological features in quantum computing provide controllability and noise error avoidance in 
the performance of logical gates. While such resilience is favored in the manipulation of quantum 
systems, it is very hard to identify topological features in nature. This paper proposes a scheme 
where holonomic quantum gates have intrinsic topological features. An ion trap is employed where 
the vibrational modes of the ions are coherently manipulated with lasers in an adiabatic cyclic 
way producing geometrical holonomic gates. A crucial ingredient of the manipulation procedures 
is squeezing of the vibrational modes, which effectively suppresses exponentially any undesired 
fluctuations of the laser amplitudes, thus making the gates resilient to control errors. 



I. INTRODUCTION 



In order to perform quantum computation (QC) effectively, it is necessary to protect our system from control 
and noise errors. While such an issue is easily dealt with in classical computation, the quantum regime involves 
intrinsic characteristics which make quantum error avoidance or error correction an important issue for QC. The 
possible errors, due to the smallness of the system under consideration, include the questionable viability of the 
^ > approximations considered and the controllability of the experimental setup used to reproduce the desired ideal 
regimes, as well as thermal noise and spontaneous emission. Two main approaches have been developed in order to 
I . overcome these problems. The use of error-correcting algorithms which theoretically enable QC to be carried out 
' even in the presence of quantum noise, and error- avoiding schemes, which produce a computational (decoherence- 
. free [||-^) subspace of the whole Hilbert space of the system, that is subject to minimal error fluctuations. Such 
schemes include topological quantum computation P-pT| , the use of cooling, adiabaticity or the Zeno effect [Q . Both 
r strategies function with resources additional to those used for ideal encoding and processing quantum algorithms. 
This extension of resources, properly treated allows the overall control or noise error to be decreased ||25[] . There is 
Oh! an analogy between the algorithmic strategies and these engineered setups. Loosely speaking, one may regard the 
error-correcting algorithm as "simulating" a corresponding error- avoiding physical process. 

In the literature there has been much interest in topological quantum computation. Ideally, it constitutes within a 
certain control procedure error-free computation. Although the proposed schemes have hitherto been hard to realize 
experimentally, the concept is worth investigating. Our particular aim here is to identify some topological features 
from the geometrical ones appearing in the ionic setup. Let us focus on the case of small errors, called here 

£. While topological QC has no dependence in e in all orders, our aim with geometrical QC is to "neutralize" our 
gates only up to a finite order in e. As an example, it is easy to visualize in our setup cancellation of the first order 
, in statistical errors in the variables the experimenter is controlling. The parameters S of the gates presented here 
■ ■ ' are engineered functions of the experimental parameters in contrast to dynamic QC, where the gate parameters are 
merely linear functions. In particular, each S can be interpreted as the area of a contour lying on a specific surface. 
This area is resilient to the first order in statistical fiuctuations of its border which the experimentalist is traversing 
The weight factor in the surface integral S, will further govern the robustness or weakness of the gates. 

~T4|,|l8|,|Tl. 



Let us initially give in an abstract way the main ideas of holonomic quantum computation (HQC) |13 
Holonomies are a generalization of the geometrical Berry phases to the case of a multiple, e.g. 2"-fold, degenerate 
Hamiltonian Hq. The quantum information is encoded in an 2"-dimensional, degenerate eigenspace C of 77o, with 
eigenvalue i^o, which is usually taken to be the lowest ground energy. The operator is considered to belong to the 
family T — {H^ — U{a) HoU^{a): a e Al} of Hamiltonians unitarily (Z^^(cr) = Z^"^((t)) equivalent to and therefore 
iso-spectral with iJg, where — for some o-g G A^- The cr's represent the classical "control" parameters used to 
manipulate the encoded states |^) G C. Let C be a loop in the control manifold M. When C is slowly traversed, then 
the evolution is adiabatic and no population is transferred between different energy levels. If \^)in G C is an initial 
state in the degenerate eigenspace, at the end of the loop it becomes |V')out = e~^^° Ta{C) \ip)in, which still belongs 
in the same degenerate subspace. The first factor is just an overall dynamical phase which is omitted in the following 
by redefining the energy levels, i.e. by taking Eq = 0. The second contribution is the holonomy Ta{C) G U{2") 
and is a result of the non-trivial topology of the bundle of eigenspaces over A4. By introducing the Wilczek-Zee 
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connection 



Al::={V\U\a)-^U{a)\v) , (1) 

where A'^'^ is the matrix element of the (Ji component of the connection, one finds TAiC) = Pexp Jq A, \24\, 

where P denotes path ordering. The set H{A) := {ryi(C);VC € A4} C t/(2") is known as the holonomy group [23|. 
fn the case where it coincides with the whole unitary group ?7(2") the connection A is called irreducible The 
transformations r^(C) for suitable Cs can be used as logical gates for HQC. 

Let us apply these ideas to the ion traps. First of all we need to identify degeneracy between two states. This is 
provided by the lowest eigenstates of a Jaynes-Cummings Hamiltonian describing two internal ionic levels "dressed" 
with the vibrational modes of the ion. Manipulating the vibrational modes in an adiabatic cyclic fashion generates a 
holonomic operator which evolves the degenerate dressed states. These control procedures, which include one and two- 
mode displacing and squeezing, can be performed by employing lasers which make coherent manipulations between 
the Fock states of the harmonic modes. Let us look more explicitly at these manipulation and the conditions that 
need to be satisfied in order to obtain degeneracy. 



II. SQUEEZING AND DISPLACING OF VIBRONIC OSCILLATING MODES 

In the following, the aim is to produce squeezing and displacing operations on the oscillating mode of a trapped 
ion. This is produced with lasers which couple the vibrational mode with two internal energy levels of the atom. In 
particular, two standing-wave lasers are to be used for squeezing. With their appropriate frequencies they will produce 
coherences between the vibrational states \n — 1) and \n + 1) which simulate a squeezing operation. In addition, a 
traveling wave will be used to insert a displacing term in the overall Hamiltonian. Let us describe briefly the setup 
first introduced in Q. 
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FIG. 1. Atomic levels with vibronic modes and laser detunings. 



A two-level trapped ion assumed to be located at a common node of two standing- wave laser fields with the frequencies 
located symmetrically about a carrier frequency ujo at which we put an additional traveling-wave field (5|. Using 
traveling waves for all lasers will produce the same result with simplifications in the experimental setup, even 
though it might make it more difficult to meet our theoretical limits 0. The resulting field <) is given by 

£;(+) (i?, t) = Ea sin(fcai?)e-'('^«-'^)*-^'^" + Eb sm{kbR)e-'^'^"+'''>'-''^'' + Eoe-''">"e-"^°'-"''° . (2) 

The Lamb-Dicke parameter rjj is defined by kjR = rjj{a + a^), where rjj = nao/Xj with ao being the amplitude of the 
ground state of the trap potential and Xj the optical wavelength. Attention is confined here to the experimentally 
viable Lamb-Dicke limit rjj <C 1 and hence the Hamiltonian can be expanded to the first order in rjj. Let us define 
a rotating frame given by the unitary transformation p' = WpU where U = exp[— i(j/a^a + ojo(Tz/2)t] and p is the 
density matrix of our system. By omitting the fast oscillating terms by assuming that v is much larger than any other 
parameters characterizing the system, the master equation for p' is given by 

where 
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H' = (T+{gaa + gho)) + -2^+ + ' W 

where a, are the annihilation and creation operators for the trap motion, (t+ :— \e){g\, cr_ := \g){e\ and cr^ := 
|e)(e| — are the Pauh matrices describing the two-level transition of frequency cjq. For the above the following 

conditions were taken; for gj = r]j^je~'^'^^ /2 with Q.j being the real Rabi frequency of the j-th laser (^o = 0) it is 
assumed that J7o «^ ^^1,2 so that terms of order 770^^0 can be neglected. 

In this Hamiltonian the excitation of the ion is coupled with two transitions in the vibrational modes. The first, 
produced by laser 1, is a lowering of the vibrational number from, for example n to n — 1, while the second, laser 
2, raises the vibrational number from n to n + 1. Coherent application of these two processes establishes Raman 
transitions between alternative vibrational states (i.e. between \n — 1) and \n + 1) ), finally giving a nonlinear 
effect. It is interesting to note that performing these transitions without the intermediate atomic levels would require 
quadratic terms in the creation and annihilation operators. Going back to a resonance frame by making the inverse 
transformation U =^ exp[~i{iya^ a + gaz/2)t] one obtains 

P = SHe)DHa)pD{a)S{e) , (5) 

satisfying dP/dt = —i[Hjc',P] with 

Hjc = lya^^a + |(7, + g{c7+a g-^^s-^)* + atcr_ e'^f"'')*) , (6) 

and e = rie'^S n = -l/2tanh"^(|5h|/|ga|), Oi = cj)^ - (t)b + 2iyt, a = -VLq / {2g)e^''\ g = Igalcoshri - |5fc|sinhri. For 
g = V, Hjc is the Hamiltonian of the Jaynes-Cummings model describing a two-level atom coupled to a harmonic 
oscillator. 



III. JAYNES-CUMMINGS MODEL AND DEGENERACY 



Holonomic quantum computation is based on cyclic adiabatic transformations which allow exchange of populations 
between states with the same energy. Degeneracy of energy thus has to be introduced in our system. The natural way 
to achieve this with the JC model is to tune the parameter g in such a way that two eigenstates of the Hamiltonian 
become degenerate. In the JC model derived in the foregoing the Hamiltonian has the ground internal state \g) 
with energy oji = —g/2 and the excited state |e) with energy L02 — g/2. The lowest energy eigenstate of Hjc is 
\g, 0) = \g) ® |0) given as the tensor product of the internal and the vibronic ground states, and it has the energy 
Eq — fiLOi. The rest of the eigenstates (dressed states) are given, for the resonant case = ti)2 — ^1 = 5, by 

|n,±) = y|(|.9,n-t-l)±|e,n)) , (7) 

with the eigenvalues En± = fi{LOi + v[n + 1) ± v^/n + 1), respectively. The resonant condition g = v also creates 
the desired degenerate condition Edeg = Eg = Eq^ and hence a two-dimensional encoding space is spanned by the 
orthogonal states 

0), |0, = y|(|5, 1) - |e, 0))}^{\^)/^ = 0, 1} , (8) 

where |0) and |1) are the qubit states. None of the rest of the dressed states can have eigenvalues the same as Edeg, 
which makes the degeneracy strictly twofold. The energy of this subsystem is Edeg = ^'^i and, in contrast to the 
rest of En±, it does not depend on the trap frequency v. The immediately higher energy state is the |1— ) state with 
energy Ei^ = n{uji + {2- V2)i^). The energy difference between it and the degenerate states gives the energy scale 
with which adiabatic changes have to be compared in order not to have mixing of our encoding states with the rest 
of the dressed states. 

The condition g — v for producing degeneracy has to be compared with the derivation of the Hamiltonian (^) 
where we omitted fast oscillating terms with frequency 2v in relation to g = u. A more plausible condition can be 
obtained if the constructed degeneracy is between the states |0, -f) and |1, — ), where the omitted terms oscillate with 
frequencies six times as large as the terms depicted in Hamiltonian (^) . In this case simulations showed that squeezing 
of the vibrational states can be faithfully produced with maximum amplitude ri « 2. Even though the results will 
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be the same as those presented here, the mathematical derivation is somewhat more compUcated, and so attention 
wiU be restricted to the simple case of the degeneracy between \g, 0) and |0, — ). In addition, having strong lasers may 
influence the internal structure of the atoms due to the Autler-Townes splitting Here we consider the case of 

weak binding limit where along the variation of the laser amplitudes the internal atomic levels employed here remain 
unaffected. 

In the previous section, it was seen that the states of the Hamiltonian H arc the squeezed and displaced states 
of Hjc- Hence, H has a degenerate subspace spanned by the rotated basis {\a,e;i) = D{a)S{e)\i) / i = 0,1}. 
This allows the degenerate states to be displaced and squeezed at will by varying the phases and amplitudes of the 
electric fields of the lasers. Note that the dressing of the internal states with the vibrational modes provides us 
with full controllability of the qubit encoded in the two-dimensional degenerate subspace simply by manipulating the 
vibrational modes without employing transitions between the internal states. 



IV. ONE-QUBIT HOLONOMIC GATES 



At this point we have all the ingredients necessary for performing holonomic quantum computation for one qubit. 
The mathematical steps now to be performed are similar to those used to calculate the holonomies for the optical 
setup With Ho being the JC Hamiltonian and M the parametric space spanned by the coordinates {x, y, ri,9i}, 
for a = X + iy and e — rie'^"^, we are able to perform iso-spectral transformations belonging to the family = 
{Her — HoW{(j) / a G M; U{(t) — D{a)S{t)} of Hamiltonians unitarily equivalent to and therefore iso-spectral 
with iJo, where Hq = H„^ for some ctq e M. With the operator A„ = U^-^U defined for U = D(a)S'(e), the 
holonomic connection of our model has the matrix elements A'^ = (ij^CTb)- Their components are given explicitly 



in the following formulas: 
|(cosh4ri - l)(2a'<a + l) - 



Aa — cosh 2ri H 
- ^sinh4ri(e^®i a^^ 



ae-*^isinh2ri -f f , = -A^, A,, = 
- e~*^^a^), eventually yielding 



Ag, = 



—iy — -^(cosh2ri — e*^i sinh2ri) 

-^(cosh 2ri — e^'^i sinh 2ri) —iy 



Ay — i{Aci — Aa) 



^(cosh2ri + e^^'^i sinh2ri) 



-^(cosh2ri + e'^i sinh2ri) 








Afl, = 



1 

i 



-(cosh4ri — 1) . 



Let us briefly mention that these connection components, A^., give rise to two particular field strength components, 
Fa^a^i = da^A^^ - da^A^^ + [A^j, A^J, analytically given by 



p 



-i 

1 



%/2e-2'-i 



for 6*1 = , 



(9) 



and 



Friy — i 



1 

1 



for 9i 



(10) 



Holonomic gates can now be constructed by traversing closed paths in A^. In particular, we shall choose two convenient 
planes in the four-dimensional control manifold and evaluate the holonomies generated by any loop lying on those 
planes. Explicitly, the loop Cj G {^iTi) ei=a gives r^(C/) = exp— io-iE/, with S/ := Jjj^^^^dxdri2e~^'^^ . The loop 



Cii e {y,ri)g^^^ gives r^(C//) = expia2T.11, with area S// := 



D{Cn 



. dydr i2e' 



-2ri 



In the above, D{Cp) with 



p = 1, 1 1 is the surface on the relevant sub- manifold (o'i,(Tj) of M whose boundary is the path Cp. These two 
unitaries, r^(C/) and r^(C//), are sufficient to produce any one qubit gate. The exponential damping feature of the 
field strength has been transported into the parameter of the holonomic gates, Sp, making them resilient to control 
errors for large values of the squeezing parameter, ri. This point will be elaborated in a following section. 
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V. TWO IONS AND TWO COLLECTIVE VIBRATIONAL MODES 



An N qubit system can be realized with a chain of N ions. The internal degrees of freedom of each ion are 
represented by a two-level system, while the external degrees of freedom are described by the collective modes of 
the crystal. These are the modes which diagonalize the potential (assumed here to be harmonic) and the Coulomb 
repulsion between the ions. For sufficiently cold crystals the ions oscillate harmonically around their equilibrium 
position. Each qubit is composed of a two-level system corresponding to the internal two-level transition of the ion, 
together with a collective mode of the motion. We produce a two-qubit gate by interactions between the collective 
modes of each qubit. This interaction is realized with a two-mode squeezing or displacing transformation. In the 
following, these transformations are constructed for two ions. 

In the case of a chain of two ions the collective modes are the "center of mass" mode, where the two ions oscillate in 
phase, and the "stretching" mode, where the two ions oscillate out of phase. Let us define qubit 1 as the combination 
of the internal levels of ion 1 with the center of mass mode, and qubit 2 as the two levels of ion 2 together with the 
stretching mode, analogously defined as in eq. (|^). Let Ri be the displacement of ion i from its equilibrium position 
and let the oscillations of the ions described by two harmonic oscillators of frequency I'l, 1^2 — \f^vi for the center of 
mass and the stretching modes, respectively. The quantized oscillations are described by the operators ai, a|, with 
i = 1, 2, which are the annihilation and creation operators of the center of mass and the stretching mode, respectively. 

Let us consider the following Hamiltonian of two ions where only the first is coupled with both collective vibronic 
modes via the appropriately detuned lasers: 

Hi = i^ialai + V2ala2 + -fal + -j-al - -[aXE^+'> {Ri,t) + a\E'^-\Ri,t)] . (11) 
The laser field acts on the first ion only, but couples it to both oscillating modes by the following configuration: 

E^+\Rx,t) = Ea sin(fcai?i)e-*('^^-'^^)*-"^» -|- E^ fiin[hRi)e-'^'^^+''^^^~">'^ , (12) 

with kjRi = rij^\ai + a\) + {a2 + a\), where index on the frequencies and on the creation and annihilation 
operators distinguishes between the two collective modes and the index on the Pauli matrices distinguishes between 
the two atoms. In the Lamb-Dicke limit we are able to expand the trigonometric functions up to the first order 
in both of the parameters yy'^-* and ry*-^-*. By moving to the rotated frame defined by p' = U^pU, where U — 
exp[—i{iyia[ai + i'2a,2a2 + "^cr] + ^a'^)t] we neglect all fast oscillating terms. These are phase factors with rotating 
frequencies i^i^2 or i^i ± The resulting master equation is given by 

^^-^[H[,p[], (13) 

with 

H[ = gPa\a2 + yf^^^I + ^-c , (14) 

where g^^^ = rf^^VLje^^'^^ /2. Acting on the second ion with the detuning of the lasers interchanged between the two 
collective frequencies, we obtain the following Hamiltonian 

H'2 = 5i''4ai + gi^^alal + h.c. . (15) 

If the two interactions described in Hi and H2 are turned on simultaneously then back to a resonant frame a two- 
mode squeezed density matrix P = AI^{()(^pi (g) p2)M{() is produced satisfying dP/dt = —i[HTot, P], where the JC 
Hamiltonian of the two atoms is given by 

Htoi = t^iajai + 1/202^2 + ^^o-z + ^^'^z + + 9(2)<^+a.2 + h.c.} , (16) 

with M(C) = exp((^a|aj — C0'ia2), C being rotated in the interaction frame defined by U and given by C = r2e**^, r2 = 

-tanh-i(|5«|/l5i'^l) = - tanh-i(|5f Vlsi'^l), 9(1) = l5i'^|coshr2 - |sinhr2, 5(2) = Iffi'^ | coshr2 - |5r^|sinhr2 
and 02 ~ —4>b + [vi + V2)i, where we have taken (pa — 0. In the same way, by interchanging the actions of the two 
laser fields on the two atoms it is possible to obtain two-mode squeezing and two-mode displacing transformations, 
the latter being defined by the unitary operator N{£) = exp(^aja2 — £,aia\). 
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VI. TWO-QUBIT HOLONOMIC GATES 



If the ground states of the two atoms \g)i and \g)2 are kept on the same energy level ujI = w^, then the two degenerate 
spaces of the two JC models, created by the conditions ^(i) = vi and (7(2) = 1^2, are also degenerate between each 
other. Hence, any adiabatic evolution will allow the mixing of these four degenerate states while they are kept isolated 
from the rest state space. The operators Ar^ = M^C)-£qM{() and Ar^ = M^C)(^N^{Ow;NiO^M{C) give rise to 

the connection components A'^J ~ (i|ylcr|:/) for \i) and |j) belonging to the tensor product basis of the two qubits 
|1) = |00), |2) = |01), |3) = |10) and |4) = |11). Hence, we obtain 
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(17) 



i(2coshV2 - 1) . 



The relevant field strength component is given by 



p = 







e*^3 









\/2 sinh 2r2 



Two-qubit holonomic gates can be produced by performing the following cyclic evolutions. By means of a loop 
Chi e (^2, J's) 02=93=0 we obtain TaCC///) = cxp-ia^^^^S/// , with E/// /^^p^^_^(ir2dr32sinh2r2. For Civ G 



(12) 

t'^'s) 92=0,03=31^/2 following holonomy is generated: Ta{Civ) — Gxp—itJi E/y , with the area given by E/y :— 
/d(Civ)^^2i^''32 sinh 2r2. In the above we have used 



-(12) 





-i 
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and 



,(12) 
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As an example, Ta{Civ) with E 



IV 



7r/4 gives the following non-trivial two-qubit gate: 
V2 







-i 
1 
^/2 



which together with the holonomies produced by loops Ci , Cii compose a universal set of transformations pE 



VII. MEASUREMENT PROCEDURE 

At the end of the holonomic procedures it should be possible to distinguish between the states |(?,0) and |0, — ) by 
a measurement. For this purpose, a 7r/2 pulse is applied to the ion: 

C/(0) = exp [ - ^ J(|e)(5|ae-** + h.c)\ . (18) 

This evolution operator will not affect the state |.g,0), but only |0,— ). In the rotated frame we are working in, the 

phase 4> is shifted to ^ = (/) — (cjq — g)t. Hence, the state jO, — ) in the laboratory frame becomes — e~*~'^*|e, 0) when 
the timing is such that — —-k 12. Direct measurement of the atomic energy levels distinguishes between the two 
logical states of the qubits. 
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VIII. CONTROL ERROR RESILIENCE 



For the one-qubit gates presented in Section IV the weight factor of the surface integrals E/ and S// is e"^''^ . Let 
us introduce errors in the displacing and squeezing parameters. For the rectangular loop of Fig. ^(a) we may take 
the borders, instead of being positioned at (x, ri), to be at (x + ei, ri + £2) for small errors ei and £2- 



X i\ 





FIG. 2. (a) The rectangular loop C on the plane of the displacing parameter x (or y) and the squeezing amplitude r\. (b) 
The cyclic evolution in the phase space (x,y): in (1) the |0) Fock state is depicted; then it is displaced along the x axis in (2), 
squeezed (3), displaced back to the origin (4) and then, to complete the cyclic evolution, it is squeezed back to the original 
state 10). 



Then the areas E/ and E// are varied by 

AE = £1(1- e-2('-i+^2)^ ^ ^g-2n 



(19) 



where its dependence with respect to both errors is given in Fig. ^. For values of ri « 2 the error £2 is exponentially 
suppressed in relation to the contribution of £1, which is linear. 

The resilience to the squeezing error is of the following nature. As has been seen, the gate parameter can be viewed 
as the flux of the field strength F (which can be regarded as a virtual magnetic field). Its amplitude decreases for 
large r\. Hence, a loop which has one boundary at a large value of r\ will include the same flux for small deformations 
of its boundary owing to the exponential dumping of the field strength amplitude. This characteristic is conceptually 
very interesting. A non-trivial bounded geometry, i.e. one whose curvature decreases exponentially from a fixed point, 
looks like topology from far away. A simple analogy can be drawn with the Aharonov Bohm effect. There, a magnetic 
field is bounded in a solenoid and a charged particle travels around it. If the particle performs a loop inside the 
solenoid, it acquires a phase given by the enclosed flux of the magnetic field. This has a one-to-one correspondence 
with the geometrical phases. But if the particle travels /ar away from the solenoid, it acquires a fixed amount of 
phase, varying only by the integer number of circulations the particle makes, which is a topological variable. 
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FIG. 3. The error in the gate parameter, AE, as a function of the displacing (ei = Aa;) and squeezing (e2 = Ari) errors. 
The latter has a much smaller (about two orders of magnitude) influence on AE. The dimensions of the rectangle are taken to 
be a; = 1 and r\ = 2. 



IX. CONCLUSIONS 

In this article, we have constructed an ionic setup where qubits are encoded in the space of degenerate dressed 
states constructed from internal energy levels of the atoms combined with their vibrational modes. The logical gates 
are represented by the holonomies acting on the degenerate states. They are produced by adiabatic cyclic evolutions 
of the oscillating motion of the ions. These holonomies are nonlinear functions of the parameters the experimenter is 
controlling. This provides a wide range of possibilities for constructing gates resilient to control errors. In particular, 
we have seen one-qubit gates which afford exponential suppression of errors. 

For two qubits the two-mode squeezings and displacings produce holonomies with hyperbolic dependence on the 
control parameters. It is possible to use two-mode control manipulations of higher order in the creation and annihila- 
tion operators than the bilinear forms we have seen with M{Q and N(£^). Even though their transformations have not 
been mathematically studied sufficiently to evaluate their holonomies, they appear promising for producing trigono- 
metric [ p9[ or even experimentally decreasing functions of the control parameters, while they are experimentally easy 
to construct pO| . 

For quantum computation the threshold of error rate before which error correction can be efficiently applied is of 
the order of 10~^. This very low limit is more likely to be achieved by advantageous manipulation strategies rather 
than by just trying by brute force to improve the controllability of our external parameters. Holonomies provide the 
advantage that the parameters of our gates are engineered functions of the experimental parameters. Furthermore, 
they are resilient to statistical errors, thus providing overall an appealing framework for quantum computation. 
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